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We suggest a way to characterize the coherence of the spht Cooper pairs emitted by a double- 
quantum-dot based Cooper pair splitter (CPS), by studying the radiative response of such a CPS 
inside a microwave cavity. The coherence of the split pairs manifests in a strongly nonmonotonic 
variation of the emitted radiation as a function of the parameters controlling the coupling of the 
CPS to the cavity. The idea to probe the coherence of the electronic states using the tools of Cavity 
Quantum Electrodynamics could be generalized to many other nanoscale circuits. 



Entanglement is now accepted as an intriguing but available resource of the quantum world. However, it is still 
unclear whether this phenomenon can survive in nanoscale electronic circuits, because an electronic fluid is charac- 
terized by a complex many-body state in general. Producing and detecting entangled electronic states are therefore 
important goals of quantum electronics. The spin-singlet (Cooper) pairing of electrons, naturally present in conven- 
tional superconductors, appears as a very attractive source of electronic entangled states. Recently, the splitting of 
Cooper pairs could be demonstrated in Y-j unctions made out of nanowires [l|-l3|. But how to distinguish between a 
singlet state broken into a product state due to decoherence and the desired coherent singlet state is not addressed 
in these experiments. Here, we show that the radiative response of such devices can reveal the presence or absence of 
entangled states. The split Cooper pairs are shown to decouple from the electromagnetic field conveyed by a photonic 
cavity (subradiance) if coherent. These findings add a new twist to quantum opto-electronics, and could be applied 
to any source of on-demand entangled electronic states. 




FIG. 1: Scheme of the CPS embedded in a photonic cavity. The CPS is made out of a double quantum dot coupled to a 
central blue superconducting contact (S) connected to ground, and two outer yellow normal metal contacts (N) biased with a 
voltage V. This circuit is placed inside a photonic cavity (represented schematically by mirrors in green). The Cooper pairs 
spread over the K and K' orbit als of the dots, sketched in blue and green respectively. The low energy level structure of the 
system allows photon emission which can be amplified through a lasing effect. 



The analogy between a beam splitter for electronic states and a beam splitter for photons calls for the use of corre- 
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lation measurements to characterize the degree of entanglement of pairs of electrons, via the noise cross-correlations 
of the electrical current (3-[iq|. The use of a double quantum dot circuit connected to two normal electrodes and 
one superconducting electrode gives a practical realization of an electronic entangler and simplifies the diagnosis of 
entanglement from transport, as originally suggested by Recher et al. [7]. Nevertheless, measuring cross-correlations 
in such a setup is a formidable task because the conditions for useful entanglement correspond to a regime where the 
system is almost isolated from the leads. In this case, the Cooper pair current is too small to yield measurable current 
fiuctuations with present amplification techniques pJJ]. These difficulties stem from the fact that probing directly a 
quantum system with transport measurements is not natural since one needs the system to be open and closed at the 
same time. On the contrary, as known in atomic physics [12|, the use of light-matter interaction is very well adapted 
to probe closed quantum systems. Here, we show how to use the coupling between electrons and photons from a 
microwave cavity, to assess the coherence of Cooper pairs emitted by a CPS implemented with a double quantum dot 
circuit. 

We consider a CPS made out of a single wall carbon nanotube in which a double dot is defined by a central 
superconducting contact connected to ground and two outer normal metal contacts biased with a voltage V. This 
CPS is inserted inside a photonic cavity (see Figure 1) which is assumed to be implemented in a coplanar waveguide 
geometry using a superconducting metal, like in the Circuit Quantum Electrodynamics architecture [13|. Very recently, 
it has been demonstrated experimentally [14,d5j that one can extend this architecture to quantum dot circuits 16Ml9|. 
bringing our proposed 
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FIG. 2: a. Spectrum of the subset S of even charge states of the CPS near 6 = 2Ar. It comprizes three spin-triplet states, and 
two states \Vi) and IV2) which are a coherent mixture of the empty state and a spin-singlet state |<S). b. Dynamics of the CPS 
near the working point S = 2Ar. We consider a bias voltage regime such that the various transitions between the CPS states 
occur together with the transfer of one electron towards the normal contacts, apart from the transition from \Vi) to the triplet 
state |T-), which occurs due to the emission of a photon towards the cavity. We have checked that the radiative transition 
from |T_) to IV2) is not active in the regime of parameters we consider. 

Our aim is to find an effect which reveals qualitatively the coherent injection of Cooper pairs. Thanks to Coulomb 
blockade, the states participating to electronic transport can be reduced to the double-dot empty state |0,0), the 
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singly occupied states |rcr, 0) = d]^^^ |0, 0) and |0, ra) = |0, 0) of the left (L) and right (R) dot, and the non-local 
doubly occupied states \r(j,r'a') = dL^^dj^^,^, |0,0), with r G {K^K^} and a G {t^i} being orbital and spin indices 
respectively (see Supplemental Material (SM), part A for details). For biasing conditions like the one depicted in Fig. 
1 with eV < A (A being the gap of the superconducting contact) Cooper pairs are injected from the superconductor 
into a singlet state |<S). The coupling to the superconductor hybridizes the \S) state with the empty state |0,0) 
forming two states which we call \ Vi) and IV2). These states can relax to a single particle state via electron tunneling 
into the normal leads, and then via a second tunneling process into the state |Vi) or IF2), which closes the usual 
operation loop of the CPS, as depicted in Fig. 2b. However, when the CPS is placed in the resonant cavity there is an 
additional transition to a triplet state (two electronic states with equal spin) through the emission of a cavity photon. 
It is the aim of this work to study these photon emission processes and show that they can be used to characterize 
the coherence of the split pairs injected into the CPS. 

The relevant energy scales of the CPS are the position e of the energy levels on each dot, which we assume symmetric 
for the sake of simpHcity, the Cooper pair coherent splitting rate teh 13, Ell, the effective spin-orbit coupling c onstan t 
Ago and the coupling between the K and K' orbitals, Ak^k' , which arise from weak disorder in the nanotube 
For teh = O5 ^so = O5 and Ak^k' = 0, the doubly occupied states cost an energy 2e = S. For teh <^ Ak^k'^ A 
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the regime S ~ 2A^, with A^ = y A^^ + ^k^k'^ most adequate for producing entangled states. It allows one 

to isolate, in the double dot even charge sector, a subset £ of five lowest energy eigenstates which are at a distance 
~ 2Ar from the other even charge states, at least. These five states include three spin-triplet states |To), and 

2Ar and the two hybridized even states [26|, [27 



|T_) with energy Etripiet 
energy Ei = \ {s - 2A^ - 
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\Vn) = yr^|0,0)+^n|5) with 



{S — 2A^)2^ for n G {1, 2} (see SM, part B.4, for the expression of the v'^s). 

Here, the definition of the spin-singlet state \S) and the spin-triplets take into account the twofold orbital degeneracy 
of each dot (see SM, part B.4). The energies of and 1^2) show an anticrossing with a width Ep = 2\piteh at 
^ — 2A^, while the energy of the triplets lies in the middle (see Figure 2. a). 

Let us first discuss the current Icps through the superconducting lead, without considering the effect of the photonic 
cavity. We call Fat the bare tunnel rate of an electron between one dot and the corresponding normal metal contact, 
while Pa denotes the probability of a state \ A) of the double dot even charge sector and Pgingie the global probability 
of having a double dot singly occupied state. In the sequential tunneling limit Fat <C ^ these probabilities follow 
a master equation -^P = M P with 
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Equation ([T]) corresponds to a bias voltage regime such that single electrons can go from the double dot to the normal 
leads but not the reverse (see SM, part C). We disregard the states |To) and |T+), which are not populated in the 
simple limit we consider. The states \ Vi) and |T_) can decay towards several different singly occupied states. The sum 
of the corresponding transition rates equals 2Vnv'^ and 2rAr respectively, which explains the presence of the factors 
2 in the 3 first columns of M. To calculate IcpSi one first has to determine the stationary value Pstat of P from 
MPstat = 0. Then, one can use Icps = R-Pstat with R = eFAr[2v^, 2^2^ 2, 1]. If the double dot is initially in the 
state \Vn)^ with n G {1,2}, there can be a transition to a singly occupied state while an electron is transferred to the 
normal leads, because \Vn) has a \S) component. Then, there can be a transition from this singly occupied state to 
\Vm), with m G {1,2}, because \Vm) has a |0, 0) component. This leads to the existence of state cycles which produce 
a flow of electrons towards the normal leads (see Figure 2.b). On the left (right) of the anticrossing point S — 2A^, 
the state |Vi(2)) is the most probable. Exactly at the anticrossing point, the states and |]^) contribute equally 
to current transport and Icps is maximum (see Fig. 3, main frame). 

Electronic spins are naturally coupled to photons thanks to the spin-orbit interaction, which exists in many types 
of conductors, including those used to demonstrate Cooper pair splitting. In our case, we will take into account a 
spin/photon coupling with the form (see SM, part F) 



a so = (a + CL^)^^\ad\^^dira = (« + CL^)hso 



(2) 



with dj^^ the creation operator for an electron with spin a in orbital r of dot i G {L, K\ and the creation operator 
for cavity photons. We use A^cr = icrA^. Inside the £ subspace, the term Rso couples to |T_) only, i.e., for 
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FIG. 3: Main frame: Current Iqpc at the input of the Cooper pair sphtter. The peak at ^ = 2Ar is due to the anticrossing 
between \Vi) and IV2). Left inset: zoom on the current peak due to lasing transitions between the \Vi) and |T_) states. Right 
inset: number of photons versus 5. We have used reaUstic parameters teh = 9 /ieV, Aso = 0.15 meV, A^/k' = 0.9 meV, 
z/o = 3.64 GHz, Q = 2i:nv^lK = 500000, D'^ = 92 ns, Vn = 72 MHz, T = 35 mK, and \\l - \r\ = 2.5 10"^ A^ ^ 0.55 MHz. 



^{1,2}, 



(T_ I hso \Vn) = iVn^^^^iXL - Xr) 



(3) 



whereas (^+[o]| ^so |^i(2)) = 0- The presence of the minus sign in Eq.(j3]) is crucial. It reveals the coherent injection 
of singlet pairs inside the CPS. If \Vn) was resulting from the hybridization of a product spin state with the |0,0) 
state, the matrix element ([3]) would depend only on or A^^. We will show below that the peculiar structure of the 
element ^ can be revealed by measuring the lasing effect associated to the transition \Vi) 

In order to study current transport and the photonic dynamics simultaneously, one can generalize the above master 
equation description by using a semi-quantum description of lasing [28]. This consists in adding inside the matrix M 
rates accounting for photonic emission and absorption processes, i.e. M ^ M + Mph with 
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and (n) the average number of photons in the cavity. The master equation must be solved self-consistently with a 
photonic balance equation 



= Pv,Wv,T. ((n) + 1) - Pv,Wt.v2 {n) 
+ Pt. [Wt_v, ({n> + 1) - Wv.T. in)] - n{{n) - {n),J 



(5) 



with = {exp{fuAJo/kBT) — 1) ^. The photonic transition rates WviT- and Wt_V2 appearing in the above 

equations can be expressed as: 
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where n is the damping rate of the resonator and Dqp the decoherence rate associated to the p < — > q resonance. 
The above description is vahd provided the transitions |Vi) |T_) and |T_) \V2) are not both resonant with the 
cavity, i.e. one does not have simultaneously S = 2Ar and 2'KhuQ = \^teh- The left inset of Figure 3 shows a zoom 
on Icps around S = Si{uo) = 2A^ — 27:hvQ + {tl^/Trhuo). A current peak occurs at (5 = Si{i'o)^ due to the lasing 
effect which involves the |Vi) |T_) transition. Note that the population inversion necessary for the lasing effect 
is achieved without any AC excitation thanks to the DC bias voltage. The tunnel transition rate from |T_) to the 
singly occupied states is larger than the transition rate from |Vi) to the singly occupied states (see Eq.([T])), which 
explains the increase in Icps while the system lases. However, for typical parameters, the current peak due to lasing 
corresponds to ~ 100 fA over a strong background of ~ 10 pA. The lasing effect is more clearly visible through the 
average number (n) of photons in the cavity, which can be measured with microwave amplification techniques [29| 
(see right inset of figure 3). 




FIG. 4: Main frame: number (n) of photons in the cavity for the same parameters as in Figure 3, as a function of the coupling 
Al, for constant values of S and Xr. The vanishing of (n) for small values of |Al — Xr\ is a direct signature of coherent injection 
of split Cooper pairs inside the CPS. Inset: maximum number of photons which can be obtained at ^ = ^i(z^o), as a function 
of the cavity frequency and \Xl — Xr\. If \Xl — Xr\ is sufficiently large, the lasing effect can be obtained for a broad range 
of Uq. 

The peculiar shape of the matrix element (j3j) is particularly interesting if one can tune independently the couplings 
^L{R)^ e.g. tune Xl while A^^ remains constant (see SM, part F.2). The main frame of Figure 4 shows (n) versus 
Xl- Strinkingly, this curve is not monotonic. When A^ is too close to A^^, the average number of photons in the 
cavity collapses, because the element (T_ | hso \ Vn) becomes too small. Two lasing thresholds appear, one for A^ < A^^ 
and one for A^ > A^^. This non- monotonic behavior is directly due to the presence of the minus signs in Eq.([3]) and 
it therefore represents a smoking gun for coherent Cooper pair injection inside the device. Note that for a typical 
CPS, the range of cavity frequency leading to the lasing effect is rather broad, because the parameter 5 can be 
tuned to Si{i'o) with the dots' gate voltages. The inset of figure 4 shows the maximum number of photons in the 
cavity, obtained at (5 = ^z(i^o), as a function of uq and \Xl — ^r\- For the parameters of Fig. 3 and in particular 
|Al — Ai^l = 2.5 10~^A^, the range 3.3 GHz < < 5.9 GHz allows one to have (n) > 40. Interestingly, for Sr > 2A^, 
the state IV2) is more probable than |Vi), but it cannot produce any lasing effect. In this case, one can imagine to 
test the existence of the minus sign in Eq.(j3]) by applying to the CPS a classical AC gate voltage with frequency 
{Efripiet — E2)/27rh instead of coupling it to the electric field of a cavity. The current through the CPS should show 
a non- monotonic dependence with respect to Xl or A^^. However, the current variations corresponding to this effect 
might be very small. The lasing effect discussed in this letter presents the advantage of providing an intrinsic and 
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powerful amplification mechanism for the coefficient (T_ | hgo This constrasts with a previous proposal where a 
cavity was not used [30| . 

Before concluding, we discuss our findings on a more general level. In principle, spin-fiip processes induced by the 
magnetic coupling Ql^r) between the cavity photons and spins in dot L{R) can also lead to the subradiance effect. 
Using a generic single orbital model for each dot, which leads to CPS doubly occupied states |cr, <j'), we find transition 
elements between the singlet state |i, t) ~ W->\) ^he triplet states It^t) ^t^id which are directly proportionnal 
to ql — QR' This suggests that the subradiance phenomenon studied here is a property of the injected states rather 
than of the spin/photon coupling mechanism. Nevertheless, in practice, gL{R) is expected to be extremely small jsij . 
One thus has to consider a spin/photon coupling term mediated by the cavity electric field, like e.g. the term of Eq. 
(|2]) caused by spin-orbit coupling. In principle, our findings can be generalized to other types of quantum dots with 
spin-orbit coupling like InAs quantum dots (see e.g. Ref. [32]), but the detailed analysis of these cases goes beyond the 
scope of the present work. In our case, we have found a coupling element ([3]) which vanishes for Ak^k' = 0. However, 
^K^K' / is not a fundamental constraint to obtain a subradiant lasing transition. Indeed, at another working 
point S ^ — 2A^, we find a third hybridized even state {Vs) = y^l — |0,0) + vs \S) which can decay radiatively to 
a triplet state \Ta). One can check that the coupling element {Ta \ hso iVs) has a subradiant form but remains finite 
for Ak^k' = 0. Here, we have chosen to discuss the lasing transition |T_) at J = Si{uo) ~ 2Ar because its 

frequency is more likely to match with the cavity frequency in practice (see SM, section D for details). 

In conclusion we have shown that the coherence of entangled states produced by a Cooper pair splitter can be 
proven by the lasing properties of the device when coupled to a microwave cavity. The idea to probe the coherence of 
electronic states using the tools of cavity QED could be generalized to many other nanoscale circuits. 

We acknowledge fruitful discussions with J. M. Raimond, J. Klinovaja and B. Braunecker. This work was financed 
by the EU-FP7 project SE2ND[271554]. 
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SUPPLEMENTAL MATERIAL 



A. Hamiltonian of the CPS 

Inside the left and right dots i G {L, i?}, an electron with spin a G {t^i} can be in the orbital r G {K, K'}^ which 
is reminiscent from the K/K' degeneracy of graphene. We use a double dot effective hamiltonian which includes the 
superconducting proximity effect due to the superconducting contact, i.e. 

^ddit = 5Z + /^sor(j)nira + Hint + ^K^K'Y^{d\j^AK' a + d\j^,JiKa) (7) 

with K = K' and K' = K. The term in teh accounts for coherent injection of singlet Cooper pairs inside the 
double dot [2o| . Taking into account the superconducting contacts with the term in teh is valid provided quasiparticle 
transport between the superconducting contact and the double dot can be disregarded (see part C of the Supplemental 
Material). For simplicity, we assume that the orbital energies in dots L and R are both equal to £, which can be 
obtained by tuning properly the dots' gate voltages. The term Hint accounts for Coulomb charging effects. We assume 
that there cannot be more than one electron in each dot, due to a strong intra-dot Coulomb charging energy. The 
constant Ago corresponds to an effective spin-orbit coupling fiij. The term Ak^k' describes a coupling between 
the K and orbitals of dot due to disorder at the level of the carbon nanotube atomic structure [22-24]. The 
hamiltonian H^^^^ must be supplemented by the normal leads hamiltonian Hieads = Xlfe^ i a^i^^^lk^a^ik^cr + h.c. and 

the tunnel coupling between the dots and normal leads Ht = Xlfe^ ri a^^lk^a'^'iTc^ + ^^^^ Qfe^cr the annihilation 
operator for an electron with spin a in orbital kr of the normal lead i G {L^R}. 



B. EXPRESSION OF THE CPS EIGENSTATES 



We now discuss the eigenstates and eigenvectors of hamiltonian ([7]) in the general case. The diagonalization of 
hamiltonian ([7|) can be performed by block in five different subpaces, the two subspaces of states occupied with 
a single spin a G {t^i}^ the two subspaces of states occupied with two equal spins a G {t^i}^ and the subspace 
comprising the empty state and states occupied with two opposite spins. 



B.l. Subspace of singly occupied states with spin a 



One can treat separately the subspace of singly occupied states with spin t and the subspace of singly occupied 
states with spin |. For a given spin direction a G {t^i}, the eigenenergies and corresponding eigenvectors of H^^^^ 
are: 



eigenenergy 


eigenvector 


£l(j — £ tee 


-l^l-a^{\Ka,0)- |0, Ka)) + ^^^'^^^ (|0, K'a) - \K'a, 0)) 


^2cr = £ -\- tee — Ar 


M - - Wl - {\Ka,0) + \0,Ka)) + {%K'o) + |i^V,0)) 


SScr = £ — tee Ar 


|.3.> - - Wl + (|i^cT,0) - %Kcj)) + ^'f'''\^^ (|0,i^V)} - |ifV,0)) 


£4a = ^ + ^ee + 


|54.) - Wl + (l^'^,0) + %Ka)) + (|0,KV)} + |XV,0)) 



with = y A^^ + ^\:/K' • ^ote that tee occurs in the expressions of the eigenenergies but not in the expressions of 
the eigenvectors because we have assumed that the left and rights dot have the same orbital energy £. 
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B.2. Subspace of states occupied with two equal spins a 



One can treat separately the subspace of states occupied with two equal spins t and the subspace of states occupied 
with two equal spins |. For a given spin direction a G {t^i}, the eigenenergies and corresponding eigenvectors of 



HliL are: 



eigenenergy 


eigenvector 
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^1-) = ^ d^''^' ^''^) - 1^*^' ^'^)) + 


\K'(T, Ka) 
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52.) = ^ {\K'a, K'a) - \Ka, Ka)) + 


\Ka, K'a) 


S-2Ar 


Ssa) — \ 




\Kcj,Kcj)-\ 




\K'a,K'(j) + ^2^^^' {\Ka,K'a) + \K'a,Ka)) 






[^ + -^) 


\Ka,Ka) + \ 




\K'a,K'a) + ^^l''' i\Ka,K'a) + \K'a,Ka)) 



with = w A^^ + (A|.y^,/2). These states correspond to generalized triplet states with spin 1. 



B.3. Subspace of states occupied with two opposite spins 



We now discuss the subspace of states comprising the empty state |0,0) and states occupied with two opposite 
spins. It is practical to first define the eigenenergies and eigenvectors of H^^^^ for teh =0: 



eigenenergy 


eigenvector 





|0,0) 


S 


1^1 ) = ^ ^' ^' t)) - \C- {K'i,Km + ^ \C- {K' I, K' t)) 


S 


l^i) = ^ ^' ^' t)) - \C+{K' i, K m + ^ \C+{K' i, K' t)) 


s 


1^2) = ^ {K K'\))- |c_ {K'UKm + ^ \c- {K if t)) 


s 


1^2) = ^ {\c+{K K' t)> - \C+{K' K t))) + ^ \C+{K i, K t)> 


8-2Ar 


2 ka) = - 1) \c-{K t, if' ;)) - + 1) \c-{K' t, if ;)> + (|c_(if t, ^ ;)) + |c-(if' t, if' ;))) 


8-2Ar 


2 1^3) = - 1) ic+(if t, if' ;)) - (t" + 1) ic+(if ' t, K ;)) + ^ (|c+(if t, if ;)) + ic+(if ' t, if' ;))) 


5^2Ar 


2 k4) = + 1) ic_ (if t, if ' ;)) + (1 - ^) ic- (if ' t, if ;)) + ^ (ic- (if t, if ;)) + ic- (if t, if ' ;))) 


5^2Ar 


2 1^4) = ( + 1) ic+(if t, if' ;)) + (1 - 1;^) ic+(if ' t, if ;)) + ^ (ic+(if t, if ;)) + ic+(if ' t, k' m 



where \C±{Tcr^T' a')) = {\tcf^t'cf') ± \t' cr' ^rcr)) / \/2. From the definition of |C±(rcr, rV)), one can view the states 
1-52), \s^) and I54) as generalized singlet states whereas 1^2), l^s) and 1^4) can be viewed as generalized triplet 
states with total spin 0. In the subspace {|0,0),|si) , \ti) , \s2) , 1^2) , l^a) , \t^) , I54) , 1^4)}, the hamiltonian ([7]) writes 
exactly: 

HliL = ^(l^i) (^il + k2) (^2| + \ti) (til + 1^2) {t2\) + - 2A.) (1^3) {U\ + 1.3) (^3|) + + 2A.) (\U) {U\ + 1.4) (^4|) 
+ V2teh i\sA) (0, 0| + |0, 0) (54I - 1^3) (0, 0| - |0, 0) (53I) 

Thus, the states I52), |ti),|t2), 1^3) and 1^4) are eigenstates of H^^^^ while the states \ss), and I54) and are 

hybridized with |0,0). This hybridization leads to three eigenstates IV2) and IV3) with energy Ei, E2 and E^. 

Figure 1 shows the energy of the different doubly occupied eigenstates (with total spin 1 or 0, i.e. given by sections 
B.2 or B.3) as a function of 5. The triplet states have energies 5^5 — 2Ar or S -\- 2Ar (blue lines), while the energies 
El, E2 and Es (red, green and pink lines) show two anticrossings at (5 = — 2A^ and S = 2A^. 



B.4. Expression of the CPS doubly occupied eigenstates for S ^ 2Ar 

In the main text, we work near the right anticrossing in Figure 1 ~ 2A^). For S ~ 2A^, a simplified expression 
of the relevant doubly occupied eigenstates of H^^^^ can be obtained by performing a diagonalization in the subpace 
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FIG. 5: Energy of the different doubly occupied eigenstates (with total spin or 1) as a function of S. 



{|0,0), \ss)^\ts)^ -Sa^y, ASa^y}. This leads to the CPS eigenstates 1^2), l^o), |^+) and |T_) discussed in the main 
text, which are defined by: 



with 



eigenenergy eigenvector 



S-2Ar 



S-2Ar 



S-2Ar 



El 



\To) = \ts) = E^{U^^- 1) \C^{KcT,K'a))\ + |C+(r t,r 



V2 



|T_) 



s/2 



|yi) = yr^|0,0)+^;i \S) 



1^2) = vr^io,o)+^^2 \s) 



and 



2t 



eh 



8tl^ + {5- 2Ar){S - 2Ar + (-l)»V8*e/^ + (<^-2A.)2) 



15) = \ss) 



rr ^ 



2^ 



(j) \C-{Ka. 



^Ei^-(n,n)) 



Note that a = ±1 stands for <j G {t^i} algebraic expressions. The state |To) corresponds to a generalized triplet 
state with zero spin, |T_) and |T+) correspond to a coherent mixture of two triplet states with finite spin, and the state 
\S) corresponds to a generalized spin-singlet state. Note that without disorder {Ak^k' = 0), \S) has components in 
|C+(r i,r t)) only and |T_) has components in \t(t^t(j) only. In the presence of disorder, \S) also includes components 
in C_(|rt, r|)), and |T_) has also components in \K(j^K'(j) and \K'(j^K(j). This enables a coupling between the 
states |'^i(2)) aiid |T_) through Hgo- 

C. BIAS VOLTAGE WINDOW 



The hamiltonian H^^J^^ can be used provided there is no quasiparticle transport between the superconducting lead 
and the dots. This requires 



- A < -eV < A 



(8) 
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In the main text, we furthermore assume that electrons can go from the double dot to the normal metal leads but 
not the reverse. This is true provided the bias voltage V belongs to a certain range which we derive below. 

(1) We recall that the double dot singly occupied states l^^^r) have energies e^^r, with i G {1, 2, 3, 4}, given in section 
B.l of the supplemental material. We assume 5 ^ 2A^, so that the states l^i) and IV2) have significant components 
in |0,0) and \S) and their energies Ei and E2 are well approximated by the expressions given in section B.4. Since 
|Vi) and 1^2) have components in they can decay towards the singly occupied state \si(j) while an electron is 
transferred to the normal metal leads (and the reverse process is forbidden) if 

-eV < ^1(2) - Si, - aksT (9) 

for i G {1, 2, 3, 4}. Above, a is a dimensionless factor of order 1 which takes into account the temperature broadening 
of the levels. 

(2) Since and|V2) have also components in |0,0), a singly occupied state can decay towards or 1^2) while 
an electron is transferred to the normal leads (and the reverse process is forbidden) if 

-eV < Si, - ^1(2) - OikBT (10) 

for i e {1,2,3,4}. 

(3) Since we assume 6 ~ 2A^, the state {Vs) has a negligible component in |0,0). Hence, it can be considered as 
a pure doubly occupied state, with energy — S -\- 2 A. The other doubly occupied states have energies S — 2A^, 
6 + 2A^, or 6 (see section B.2 and B.3). The doubly occupied states (including | V3)) can relax to the singly occupied 
state \sicr) while an electron is transferred to the normal leads (and the reverse processes are forbidden) if 

-eV <5±2Ar- €ia - aksT (11) 

and 

-eV < S - - aksT (12) 

Since E2 < S - 2A^ < Ei,S < S ^ 2A^ < Es and Si^ < e2a,£3a < ^4^, the combination of Eqs. ([8]), (j9j), ([lO]), ([TT]) 
and (p!2]) yields the constraint 

-A < -eV < min(£;2 - ^4cr, ^la - ^1) - aksT 
We note 6 = 2s. For S ~ 2A^, one has 

E2 - e^a = -2Ar - tee - 8tl^ ^ {S - 2Ar)^ = Si, - Ei - 2A, 

We conclude that we have to satisfy 

- A < -eF < -2 A, - tee - \ ^8tl^^{S-2ArY - aksT 

with A the BCS gap of the superconducting contact. Since the temperature T is much smaller than A^ in a typical 
experiment, we simplify this criterion as 

-A < -eV < -2Ar - tee - ^ ^St^^ + ( J - 2 A^)^ 

With the parameters of Fig.2, tee <C A^, (5 - 2A^ and A the BCS gap of NbN, this gives 1.8 meV < eV < 3 meV. 

To populate the state iVs) and the triplet states other than |T_), one needs to have a transition from a singly 
occupied state to one of these states. In the regime we consider above, this is not possible since this would require an 
electron to go from the normal metal leads to the double dot. Therefore, the state \Vs) and the triplet states other 
than |T_) are not active. In contrast, the state |T_) can be populated due to lasing transitions \T-). 
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D. EFFECT OF THE SPIN-ORBIT COUPLING 



In order to discuss the effect of the term hgo appearing in Eq. (2) of the main text, it is practical to redefine the 



eigenvectors of H^^^^ in the subspace of states occupied with two equal spins as 





eigenvectors 


s 


\rr \ — '^it)- 


S2i) a- 'S'2t) — 









\rr. \ — ^^it)- 










s 






s 


\T2-) = ^ 




6-2Ar 


T+) already defined in section B.4 


S-2Ar 


T_) already defined in section B.4 






6 + 2Ar 





with 



an 



Al — Ai? =F -^(Al + Ai?) 



For S ^ 2Ar and the bias voltage conditions considered in section C, the states l^i) and 1^2) can be populated but 
not IV3). Only two of the triplet states, namely |T_) and |Tj,), are coupled to l^i) and IV2) by hgo- One has 



{T_ I hso \Vi(2)) = i^i(2) (Al - Xr) 



which corresponds to Eq. (3) of the main text, and 



{n\hso\v. 



(2)/ 



i^i(2)AL(A^-Ai) 



A. AL(Ai + X%) + -^{Xl - Xnr 



(13) 



(14) 



The couplings to |T_) and |T^) are both subradiant since they vanish for Al = Xr. For Ak^k' = 0, the coupling 
between |T_) and |^i(2)) vanishes, whereas the coupling between |Tj,) and |^i(2)) persists: 



hm ^{ib\riso\Vu2)) = > ^ — 



(15) 



We conclude that it is not necessary to use a finite Ak^k' to obtain a coupling between |Vi(2)) and a triplet state. 
It is nevertheless impossible to use the transitions |Vi(2)) |Tj,) for lasing since |Tj,) is higher in energy than l^i) 
and \V2). 

For S ^ — 2A^, the state IV3) can be populated. When Ak^k' = 0, we find that only one of the triplet states, 
namely |Ta), is coupled to IV3), with a coupling element 



lim {Ta\hso\Vs) 



When Ak^k' 7^ 0, one has: 



{Ta\hso\V3) 



A. AL(Ai + A2,) + ^(Ai-A^)2 



(16) 



(17) 
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These equations are analogue to Eqs. ([T5]) and ([T4|) . In principle, if the cavity frequency is matching with Es — 
6 ~ 2Ar + V^teh^ there can be lasing transitions from IV3) to \Ta) since iVs) is higher in energy than \Ta). The 
\Vs) ^ \Ta) transitions are subradiant since (T^l hgo \ Vs) cancels for = A^^. Importantly, (T^l hgo \ Vs) remains finite 
for Ak^k' = 0. Therefore, Ak^k' 7^ is not a fundamental constraint to obtain a subradiant lasing transition 
in our system. Nevertheless, in practice, the frequency of the \ys) ^ \Ta) transition is not likely to match with the 
cavity frequency because A^/k' 1 ^so ^ ^ttuq is expected [111, We have chosen to discuss the lasing transition 
^ |T_) at 5 ~ 2Ar because it corresponds to a frequency of the order of \^tehi which is expected to be much 
smaller [Jj]. 

E. OTHER POSSIBLE LASING TRANSITIONS 

• In principle, there can be radiative transitions from |T_) to 1^2), which are taken into account by Eq.(7) of 
the main text. However, the lasing threshold corresponding to this transition is not reached in the regime of 
parameters we consider, because the population of state |T_) is negligible for S > 2A^. 

• Since the parameter S depends on the dots' gate voltages, cavity photons can also couple to the CPS through the 
operator Sdiff = Er,a,r',a' \r(yy(y') {tg,t'g'\. One finds {Vi\5diff \V2) = V2teh/ V^tl^^^{8^2A^ . Hence, 
in principle, there can be lasing between the states |Vi) and IF2). Nevertheless, this can be avoided by using 

= 2TihvQ < 2\/2teh or by tuning properly J, so that Eq ^ Ei — E2. 

• In principle, spin-orbit interaction can also induce lasing transitions inside the CPS singly occupied charge sector, 
corresponding to energy differences 2A^, 2A^ + 2tee, 2A^ — 2tee and 2tee- Since the scale A^ is expected to be 
much larger than £^0, only the transitions with frequency tee/Trh can possibly match with the cavity frequency 
uq. However, it is rather unlikely to have such a matching in practice. In the main text we assume uq 7^ tee/Trh. 
This criterion can be checked experimentally by extracting tee from the data. 

F. EVALUATION OF THE SPIN/PHOTON COUPLING IN A CARBON-NANOTUBE BASED 

QUANTUM DOT 

In this section, we estimate the spin/photon coupling Xl{r) which can be obtained in the single- wall carbon-nanotube 
based quantum dot L{R) thanks to spin-orbit coupling. 

F.l. Electronic wavefunction in the absence of inter-subband coupling elements 

The position it of an electron on the nanotube is marked with a longitudinal coordinate ^ and an azimuthal angle 
i.e. it = + Rcos[(p]lt + Rsm[(f]'^ with R the nanotube radius. We write electronic wavevectors under the 
form 

1^) = e^"^ \^) \a) (18) 

where k is the electronic circumferential wavevector, \a) denotes the spin part of the wavefunction and is the ^- 
dependent orbital part, which has a structure in sublattice space. We use the spin index a G {t^i} the sublattice 
index r G {K^K'} or equivalently a G {+, — } and r G {+, — } in algebraic expressions. For a zigzag nanotube, \(p) 
is an eigenvector of 

d 

HsWNT = hv(TKSi - ^-52^) + (JtAI^Sq + (TtAI^Si - AgSi + V(0<5o 

with y(^) a longitudinal confinement potential and jso, si, 52, 53) the identity and Pauli operators in sublattice space. 
We have used the same conventions as in Refs. [22, 39j to write Hswnt- The motion of electrons along the nanotube 
circumference is quantized, i.e. 
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with the subband index and rj G { — 1,0,1} a parameter which depends on the nanotube chiral vector. We have 
introduced in the above hamiltonian intra-subband spin-orbit couphng terms in A^^, A^^ and a spin-independent term 
in A^, which are derived e.g. in Refs. [ssl, [s^. The constants Aj^ and A^^ are first order in the atomic spin-orbit 
interaction Vso and the nanotube curvature R~^^ whereas A^ is proportional to R~^. 



For a problem uniform in the ^ direction {V{^) = 0), the eigenstates = ^rNkah)^^ above hamiltonian 



satisfy Hswnt V^?, 



N,k,a 



1 



, with k the electrons longitudinal wavevector. One can check: 



hvFTK-\-aTA'^^ — Ag—ihvFk 



and 



Er,N,k,a = CrrA^^ + bJ {hVpTK. + (7t/\\^ - /^gf + {hvphf 



(19) 



with 6 = ±1 for the conduction/ valence band. In the following, we use 6 = 1. In order to define a quantum dot, we 
take into account a rectangular confinement potential 

Vconf for ^ < 
V{i) = { for < ^ < L 
Vconf for ^ > L 

We obtain confined electronic states |?/^) = |^r,Ar,n,cr), with n the index corresponding to a longitudinal confinement 
of electrons. More precisely, one has 



N^^k^^ r,N,-i~k^,a for ^ < 



V2^ 



1 



C^^^^k-^ Ur,N,k-,a D^^^-ik^^ Ur,N,-k-,a foTO<^<L 

V2^ \ 1 / I 1 / 



1 



with ki>0 and Re[^f ] > 0. The wavevectors ki and ki can be obtained from the energy-conservation condition 

Er,N,k^,a = ^r,iV + ^conf (20) 

the constraint 



and the secular condition 



exp {i2k'lL) 



n- <k^ < (n + 1)- 

1j Lj 



N,-ik'l,a 



r,N,k^,o^ 



which results from the continuity of \^|^T,N,n,a) at ^ = and ^ = L. The constants Ar^N,n,a^ Br^N,n,a^ 
Cr,N,n,a^ and Dr^N,n,a cau be obtained from the continuity of \^pT,N,n,a) and its normalization condition, i.e. 



F.2. Effect of the inter-subband coupling elements 



In this section, we discuss the coupling between electronic spins and cavity photons, mediated by the electromagnetic 
field of the cavity. 
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We assume that the nanotube is parahel to the cavity central conductor. We take into account the interaction 
of electrons with the vector potential ^ of the cavity treated in the Coulomb gauge fssl, |4o| . We quantize the 
field ~1 in terms of the photonic operators [4l[ |42[. This gives a coupling operator 



We have used above sin((/:)) = — /i_)/2z, where the operator /i± = e^*'^ increases/decreases the index N. 
The subband and spin subspaces are coupled by a term [33|, [3^ 

Winter = -^L^2^(cr-/i+ - Cr+/i_) 

which is first order in the atomic spin-orbit interaction and the nanotube curvature (spin-orbit inter- 
action in carbon nanotubes was also discussed in Refs. I35l439[ ). Here, a± is the operator increasing/decreasing 
the spin index in \a). 

Note that the terms H(^^^^ and i^f^^er ^PP^Y f^^^ wavefunctions [see Eq. (p!8|) ] 

|^r,iV,n,a) = 6*^^+^)^ |V^r,iV,n,a) ^ W) 

One can perform a first order perturbation of these wavefunctions by H-^^^^: 

n' 

with 

^ (^r,A^-a,n^a| -^fnter |^r,Ar,n -a) 
'^r],r,N,n,n' 771 771 

J^T,N,n,-(7 ~ J^T,N-(7,n' ,a 



^r,N,n,a ) for =^ and a This yields intra-subband 



The term i^-^^er couples the perturbed wavefunctions 

spin/photon coupling elements, which write for the lower subband = which we consider: 

Ar = {^^r,iV=0,n,-l(^) ^i^ter ^r,Ar=0,n, + l 

-pfULfx- I A+ *K 1 Ta- - A+ *^ 

with P = iehVrms / ^^^ef f Ri^od- 

One can check that Xk = = — ^A is purely imaginary due to the assumptions used above [zigzag nanotube 
and — ^) = ~ One thus obtains the intra-subband spin/photon coupling term of Eq.(3) of the 
main text. For L = 100 nm, R = 1 nm, r] = d = 5 /im, Kms = 4 /iV, = 3.64 GHz, Econf = hvp/^R, 
rrieff = h\n\/vF = 4.4 10~^^ kg, and using the parameters A]^ = —0.08 meV, A^^ = —0.32 meV, and 
5.7 meV taken from Ref. [sl we obtain A 0.4 MHz. 

In order to obtain a tunable spin/photon coupling, one can insert in one of the dots gate voltage supply a 
tunable capacitance made out of a single electron transistor (SET). The electric field seen by the nanotube can 
be modulated electrostatically by placing the SET in the blockaded or transporting regimes. This allows one to 
vary the couplings Ax,(^) for dot L{R). 

In the case XiK 7^ A^^', with i the dot index added in the main text, the results of the main text can be 
generalized straightforwardly by using 

(T_ I hso \Vn) = ^n ^^^^ {>^LK + XlK' - Xrk - XrK') (21) 



